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27ms paper is concerned with a study of traffic and message delay in 
a ring switched data transmission system. The system, by asynchronous 
multiplexing and data storage, shares transmission facilities among' many 
users. It is the random component of message delay due to buffering that 
is the focal point of our study. 

The basic configuration of the system is a ring connecting stations 
where traffic enters or leaves the system. A mathematical model of the ring 
is developed which accommodates an arbitrary number of stations and 
any given pattern of traffic between stations. Studied in detail is the uniform 
traffic pattern in which each user is identical and communicates equally 
to all others. Intrinsic to the model is a recognition of the bursty nature 
of data sources. Other factors that are taken into account are line and 
source rates as well as the blocking of data into fixed size packets. Formulas 
are derived from which average message delay induced by traffic in the 
ring can be calculated. 

The results of the study are presented in a set of curves where normalized 
delay due to traffic within specific system configurations is plotted as a 
function of the number of stations and source activity. The delay here is 
normalized to average message lengths. An important parameter of these 
curves is the ratio of source rate to line rate. The results show that, in 
certain quite reasonable circumstances, the delay is less than two average 
message lengths. 

Rings of 10, 50, and 100 users have been simulated on a digital com- 
puter. Data obtained from these investigations are presented and compared 
to the theoretical estimates for line busy and idle periods and message 
delay. The results of simulated average message delay show that, for inter- 
station link utilizations of about GO percent, the difference between the 
theoretical estimates and experimental observations is small. 
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I. INTRODUCTION AND BACKGROUND 

In a recent paper 1 J. R. Pierce has proposed a data communication 
network in which users are connected in a ring or loop topology. In 
this paper we study the behavior of this network. We examine the 
relationship between source and line utilizations and the message 
delay within the system. We propose mathematical models of station 
behavior and, by making use of reasonable input data traffic models, 
predict the average delay as a function of network parameters. A 
principal result of our study shows that, in many cases of interest, 
average delay in storage is less than two message lengths. 

The network is buffered, operates on a distributed control philosophy, 
and user entry is gained by asynchronous multiplexing into the line 
bit stream. Other systems have used the ring topology. For example, 
IBM offers a synchronous, nonbuffered system 2 in which a central 
controller monitors the loop and allows access by the users. Buffered, 
centrally controlled systems have been proposed by W. D. Farmer 
and E. E. Newhall 3 and A. G. Eraser."' In the first a computer controls 
several peripheral devices, while in the second, several computers are 
interconnected. 

The characteristics of data sources and the requirements of users 
make the technique of message switching applicable to data com- 
munications. One important characteristic is that data sources are 
often bursty; i.e., relatively short sequences of bits followed by long 
pauses. A basic need of many data customers is rapid response to data 
bursts. One way to meet this need is by devoting a line to a source for 
a long period of time. However, because of the long pauses between 
data bursts the line will be underutilized. On the other hand, dropping 
the line between data bursts may be inefficient because of long setup 
times. For the line switching techniques that are currently used, the 
setup time will often be longer than the holding time. 

Message switching provides a way of efficiently using the line and 
obtaining rapid response. The idea is to asynchronously multiplex 
several sources onto the same line. For example, in the present system, 
data from each source are formed into fixed size packets and supplied 
with a header. The header contains source and destination addresses 
as well as bookkeeping information. Packets from all sources are buffered 
and fed onto the line according to some scheduling algorithm. (One 
such algorithm will be seen presently when the detailed operation of 
the ring is described.) Clearly as more packets seek access to the same 
line, more storage is required. The storage requirements and the at- 
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tendant buffering delay are important characteristics of system opera- 
tion. It is precisely these characteristics that are the focal point of our 
study. 

A sketch of the basic ring system is shown in Fig. 1. As indicated, 
traffic flows in one direction around the ring from station to station. 
The stations are indicated as B-boxes in Fig. 1. For purposes of ex- 
planation let us begin with the operation of the B-box before con- 
sidering the other components of the loop (A-box and C-box). The 
data source is connected to the B-box where its output is formed into 
fixed size blocks or packets and supplied with header. A message from 
a source may consist of several packets. If the line is free a packet 
is multiplexed on the line immediately. The packet is then passed 
from B-box to B-box to its destination. At each B-box on its itinerary 
the address of a packet is examined to determine whether the packet's 
destination is at that particular B-box. This examination entails a 
fixed delay for each B-box which can be calculated given the source 
and destination. It is, however, the random delay encountered by the 
last bit of a message before it gets on the loop that commands our 
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Fig. 1 — Pierce ring. 
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attention. A fundamental property of the system is that traffic presently 
in transmission has priority over traffic seeking entrance to the ring. 
If the line at a B-station is busy with information packets passing 
through it, the packets produced at that station are buffered until 
the line is free. The reading onto the line of a message consisting of 
more than one packet will be interrupted when packets from another 
station pass through. The reading of the message is resumed from the 
point of interruption when the line is free again. 

In order to explain the mechanism of multiplexing packets on and 
off the line, it is helpful to draw an analogy between the ring and a 
conveyer belt (see Fig. 2 for an illustration). Time slots into which 
packets may be placed circulate around the loop. The A-box insures 
that synchronism is maintained (see below). At the beginning of each 
time slot is a marker indicating whether the ensuring packet slot is 
empty or full. The B-box senses this marker and acts accordingly. 
In a full packet, address bits follow the occupancy marker. If a B-box 
senses its own address the packet is removed from the line and sent 
to its destination. This same B-box may take advantage of the empty 
slot to feed its own packet on the line. Of course if a packet slot is full 
and destined for a source beyond a particular B-box, then the line is 
momentarily blocked for that B-box. 

While there are many B-boxes in a ring, there is only one A-box. 
The A-box has two basic functions. The first, as mentioned earlier, 




X= FILLED PACKET 

Fig. 2 — Conveyor belt. Packets enter at 1 and leave at 2. 
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is synchronization of the ring. The second function is that of preventing 
the buildup of traffic in the ring due to undeliverable packets. The 
header of each packet passing through an A-box is marked. If a packet 
tries to pass through an A-box a second time it is either destroyed, 
creating an empty packet slot, or sent back to its destination. Sending 
a packet back to its destination is done simply by interchanging source 
and destination addresses. In this way a busy signal is provided. 

The C-box shown in Fig. 1 provides interconnection of rings (see 
Fig. 3). Packets destined for a station outside a particular ring have 
addresses indicating this and are picked off by the C-box in exactly 
the same way that intraring traffic is picked off by B-boxes. This traffic 
is buffered and multiplexed onto the next ring in the same way that 
traffic from a local station is multiplexed on a ring. Since traffic already 
on a loop has priority, inter-ring traffic will suffer some delay. A likely 
realization of the C-box suggested by W. J. Kropfl 5 is the tandem 
connection of buffer and B-boxes shown in Fig. 4. 
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Fig. 3 — Three-stage hierarchy of Pierce ring. 
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In the remainder of this paper a mathematical model of the ring 
will be developed and analyzed. The results of this study will be pre- 
sented in the form of sets of curves which illustrate the behavior of 
normalized delay as a function of either the number of stations in the 
ring or the utilization of the source. 

In order to carry this analysis forward it was necessary to make 
several assumptions and approximations. Thus in order to verify our 
results and to refine our model, a simulation program was developed. 
Simulation results have been obtained for rings ranging from 10 stations 
to 100 stations. These results are compared to the results of analysis. 

II. GLOSSARY OF TERMS 

C b — Line rate in bits per second. 

C p — Line rate in packets per second. 

b t — Bit rate during the activity period of the ith source. 

B p — Number of information bits per packet. 

H — Number of header bits per packet. 
1/X,— Average duration of idle period of ith source in seconds. 
1/ni — Average duration of active period of ith source in seconds. 

Q { — Average number of packets per message at ith source. 

ji — Ratio of source packet rate to the line packet rate. 

Ui — Utilization of source i [see equation (3)]. 

0, — Intensity of source i [see equation (4)]. 

r, — Average number of packets per second from source i. 

N — Number of stations in ring. 
P u — Portion of traffic from station i to station j. 

R%— Traffic passing through station k in packets per second. 

R k — Traffic out of station k in packets per second. 

U* k — Line utilization as seen by station k [see equation (9a)]. 

Q* k — Line intensity as seen by station k [see equation (10)]. 

U k — Line utilization after station k [see equation (9b)]. 

G fc — Line intensity after station k [see equation (10b)]. 
1/At — Average duration of line idle period in seconds as seen by 

station k. 
l/M% — Average duration of line busy period in seconds as seen by 
station k. 

III. SYSTEM MODEL 

3.1 Source Input Model 

In this section a mathematical model for a ring is presented thereby 
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Fig. 4 — Tandem connection of B-boxes for ring switching. 



laying a foundation for succeeding sections where the model is analyzed. 
We are primarily interested in the traffic characteristics of the ring 
and the attendant delay. In order to focus on this aspect of ring opera- 
tion, we assume that there are no equipment failures or transmission 
errors so that, once on the ring, a packet is ultimately delivered to its 
destination. Thus the action of the A-box in destroying packets or 
providing a busy signal is not part of the model. At present the study 
is confined to an analysis of a single ring. Finally the study was predi- 
cated on light to moderate loading of the system. In a heavily loaded 
system long queues of messages form and the response time of the 
system may be excessive for data applications. 

Throughout the analysis the parameter N designates the number 
of B-boxes or stations that are on the ring. The capacity of the line 
connecting the stations to one another is designated as C b bits/second. 

To each B-box is connected a source, which as we have noted above, 
is bursty in nature. The output of the z'th data source is modeled as 
consisting of alternate idle and active periods. During the latter, 
transmission is at a constant rate of 6, bits/second. The durations of 
the source activity and idle periods are assumed to be exponentially 
distributed and statistically independent of one another. Successive 
idle and busy periods are also independent of one another. The average 
durations in seconds of the active and idle periods of the source connected 
to the z'th station in the loop are denoted as l//z, and 1/X, respectively. 
Studies of computer user statistics 6,7 indicate that the foregoing is a 
reasonable, if somewhat simplified, model of a data source.* 



+ An initial study of this model of the data stream is due to R. J. Pile. 8 
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Each message (bits generated in a source activity period) is bundled 
into an integral number of fixed size packets. The maximum number 
of information bits in each packet is B„ . The number of header bits 
that accompany each packet is denoted by the parameter H. In 
Appendix A it is shown that the number of packets in a message is 
geometrically distributed with mean 

Q<=z } B /h s , i=l,2,---,N. (1) 

1 - exp (—Bjjii/bi) 

During an active period of a source, the rate at which packets are 
produced is /*,(?, . The rate at which packets can be transmitted on 
the line is C b /(B P + H). An important quantity in our consideration 
is the ratio of source packet rate to line packet rate, 

7,=^, i=l,2,---,N. (2a) 

It may well be that in some applications a typical message consists 
of many packets, i.e., b ,■//*,- » B p . In this case we have 

7<£*^- (2b) 

Other source parameters that are important in our study can be 
derived. The utilization of source i or the fraction of time that source i 
is active is 

u^r—f, i = l,2, .-.,#. (3) 

X, + M. 
The intensity of source i is defined as 

0,=^, i=l,2, •••,#. (4) 

The average number of packets/second transmitted from source i is 

>•■ = iA XT/ = M.Q.u. • (5) 



IV. LINE TRAFFIC 

4.1 Line Utilization 

From the way in which traffic is multiplexed onto the line it is clear 
that message delay is dependent upon the line traffic. As a prelude to 
the calculation of delay, the relevant characteristics of line traffic are 
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considered in this section. A precise mathematical characterization 
of line traffic is extremely difficult. In fact, for reasons that will be 
explained presently, we encounter the most difficulty in this phase 
of the analysis. 

We begin by calculating the average line utilization at each point 
in the ring. The basic assumption in this calculation is that a con- 
servation law holds so that, over a sufficiently long time period, the 
average packet rate into a station is equal to the average packet rate 
out of a station. The traffic into and out of a station includes data to 
and from a customer connected to the station as well as line traffic. 
The implication is that there is no continuous buildup of packets in 
storage at a station, which is as it should be for normal ring operation. 

Let the number of packets/second emanating from the source con- 
nected to station i (i = 1, 2, • • • , N) be designated as r, . In terms of 
parameters defined earlier /-,- = Q,[l/K, + 1/m.T'- Pa is defined as 
the portion of traffic originating at station i that is destined for station 
j with P it = 0. The average number of packets per second going from 
station i to station j is Pur, . All of these packets pass through each 
station on the ring between stations i and j. The average number of 
packets per second from station i passing through station k is given by 



-, £ Pa + '•■ £ Pa if 1 < i < fe > * * N 



R% = 



-, £ Pa 

\ t Pa 



-i £ Pa 



otherwise. 



if 1 < i < k, k = N 



if i = 1, k 5* N 



if k + 1 < i ^ N 



(6) 



The total volume of traffic passing through station k is given by 

Rt = £ R% . (7) 

r = l 

The total volume of traffic out of station k, including traffic from the 
local source, is 

R k = R* + r„ . k = \,2, ■■■ ,N. (8) 

Perhaps the distinction between R k and R* k here can be emphasized 
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by referring to Fig. 5. Here a B-box is shown as being split into its 
two functions, viz., taking data off the line and reading data onto the 
line. At point Z the line carries all of the traffic out of station k at an 
average rate of R k packets/second. At point Y the average traffic rate 
seen by the local source as it attempts to multiplex data on the line 
is R* k packets/second. (Throughout the remainder of the analysis an 
asterisk on a line traffic parameter denotes a quantity as seen by the 
location station after message deletion.) Since the packet capacity 
of the line is C„ , the line utilization as seen by the source at station k is 

U% = R*JC V . (9a) 

The line intensity at this point is 

6? = RV(C P - Bt). (10a) 

The utilization and intensity on the line after station k are given by 

U k = R k /C p (9b) 

G, = R k /(C p -R k ). (10b) 

In the case where ring traffic is symmetric, i.e., 

P, = 1° i = j (ID 

[l/(N - 1) i* j 

and r t = r,i = 1,2, ■■■, N, these expressions simplify greatly. We have 

f N 
2 



Rf =r 



- 1 



(12) 



4.2 Busy Period 

We turn now to the calculation of the average duration of idle and 
busy periods on the line. An exact calculation of either quantity is 
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Fig. 5 — Station decomposition. 
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difficult for all but very specific cases. The reason for this difficulty 
can be seen by examining the mechanism for multiplexing traffic on 
and off the line. Packets destined for a particular station may be inter- 
spersed in a sequence of contiguous packets. This sequence is broken 
up at random when the packets are delivered. On the other hand, 
the packets with the same destination may be concentrated at the 
beginning of a long sequence. In this case the long sequence will be 
scarcely affected when the packets in question are taken off the line. 

In the remainder of this section are presented two approaches to 
the problem of calculating the average duration of the line busy period. 
While neither calculation is exact, both take into account the various 
factors involved, and, in those special cases where the answer is known, 
the same correct result is obtained. The basic difference between the 
two methods lies in the kind of line traffic that they are designed to 
model. As we proceed in the calculation this will be pointed out. 

Both of these approaches draw on a result from storage theory 
concerning the data sequence out of a buffer. Suppose that n sources 
feed into a common buffer. The sources turn on and off at random 
transmitting at a constant rate during the active period, and data are 
read out of the buffer at a constant rate. The output will consist of 
successive idle and busy periods. Suppose that the durations of input 
idle periods are exponentially distributed with mean 1/a, , i = 
1, 2, • • • ,n. By definition, the probability of an input idle period 
terminating in any incremental interval dl is a, dl. An output idle 
period terminates when any one of the input idle periods terminate. 
Thus the probability of an output idle period terminating in any 
incremental interval dl is ^" =] a, dl; consequently the duration of 
the output idle period is exponentially distributed with mean 1/^"=, a,- . 

In both calculations the line traffic intensity [see equations (10a) 
and (10b)] is assumed to be equal to the ratio of the average durations 
of the line busy period to the line period. It is implicit in this assumption 
that the active and idle periods on the line are independent of one 
another and that they are stationary in the mean. 

The first method that calculates the average duration of the busy 
period looks at each station individually. In equation (6) an expression 
for R% , the number of packets per second from station i passing through 
station k, is given. If this were the only source active, the line intensity 
at station k would be 

ti* 
B " " C - Bit 
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The average length of a message from source i is Q, packets. We take 
the average duration of the idle period of the packet stream in seconds 
from the 2th source as 



1/Aft = Q t /C 9 Q% = M ^ m 



(13) 



Now the assumption is made that the durations of the idle periods 
in the packet streams are exponentially distributed. To find the average 
duration of the idle period of the packet stream as seen by the source 
connected to station k, we call upon the result from storage theory 
quoted earlier. The duration of the idle period is exponentially 
distributed with mean 



1/Af = 1 / £ Aft 



(14) 



The average duration of the busy period of the line in seconds as seen 
by the source at station k is given by 



1 
Mi 



A? 



Combining (13), (14), and (15) we obtain 



_1 
Ml 



y C v Rjk 
h Q,(C, - Bit) 



(15) 



(16) 



When the ring traffic is symmetric [see equation (11)], these results 
simplify. We have 



A * fy 1 - 71,(0 - 1)/(N - 1)) ' 



(17) 



where n = n { , = 0, , and u = w, , i = 1, 2, • • • , N. In many cases 
of interest y is small so that the term yu((j — l)/(N — 1)) in the 
denominator of (17) contributes little and may be ignored. In this 
case we may make a convenient approximation: 



A A * ^ fiu 



(f 



and 



k = 1, 2, 



N 



(18) 



1 

Mt 



1 



b-4 



- yu[f - 1 



(19a) 
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This expression may be put in the alternate form 

l/M* = 7 (i + e*)/ M . (i9b) 

The key assumption in this calculation suggests the region of ap- 
plication of this model. In the packet stream from station i passing 
through station k, at average rate R,% , it is assumed that the messages 
stay together. Thus the average busy period of this stream is Q, packets 
long. The model is not valid when the messages from individual sources 
are broken up in the act of multiplexing. Thus the model should hold 
best when the sources emit short messages or when line utilization 
is low. 

The foregoing method of calculating the duration (length) of the 
average busy period on the line essentially ignores the ring and treats 
each station as a separate entity. In contrast the second method makes 
explicit use of the ring structure. The algorithm begins with the as- 
sumption that the length of the line busy period at the input to a 
station is known. Based on this assumption the length of the busy 
period on the line out of the station is calculated. The process continues 
all the way around the ring until one returns to the starting point. 
The ring is closed by setting the duration of the final busy period equal 
to that of the initial busy period. 

The algorithm for calculating the change in the busy period is best 
explained by referring to Fig. 5. The line intensities at points X and Y 
are Q k - t and B% respectively. The assumption, fundamental to this 
approach, is that whole busy periods are deleted from the data stream. 
While there may be fewer busy periods at point Y, the average length 
of a busy period is the same from X to Y. The average durations of 
the idle periods are related by 

_L = §*zi_J_. ( 20) 

a? e A * a,_, ^ J 

The durations of the idle periods between Y and Z can be related by 
calling upon the previously quoted result from storage theory. If the 
length of the idle period at point Y is assumed to be exponentially 
distributed with mean 1/A* , and the length of the source idle period 
is exponentially distributed with mean l/\ k , then the duration of the 
idle period at point Z is exponentially distributed with parameter 

A, = A* + X, . (21) 

The average duration of the busy period at point A' is given by 

k - 1- ^ 
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One can continue in this fashion until the starting point is reached. 
In the case of a ring with a symmetric distribution of traffic the 
solution simplifies considerably since, by assumption, we have 

1 J_ 

#*-i ~ M k ' 



Substituting into equations (20), (21), and (22) we have that 

i e - e* = e - e* 

M k ~ X iid ' 



(23a) 



where \ = X,- , and 6 = 0, , i = 1, 2, • • • , N. This result can be put 
into the form 



1 7 



M k (X + |.)[1 - yu(N/2)][l - yu(N/2 - 1)] 



(23b) 



The key assumption in this calculation is that entire busy periods 
are destined for a single station. As the line traffic begins to build up, 
messages from different stations will tend to cluster in the same sequence 
and the validity of the model is weakened. The model should hold 
well when messages are multiplexed into light to moderate traffic. 

V. DELAY CALCULATIONS 

In this section we shall consider models for calculating message 
delay. As mentioned earlier we are primarily interested in calculation 
of message delay that is induced by traffic in the ring. Other delays 
such as propagation and processing delays are invariant with traffic 
intensity and are fixed for a given implementation of the ring. 

Two separate approaches to the calculation of message delay have 
been considered. These approaches are complementary in the sense 
that they apply to different kinds of source traffic while neither approach 
is applicable to the whole range of source traffic. The exponential 
on-off model of the source as presented in the foregoing is difficult to 
handle analytically. The difference in the two approaches lies in the 
way this exponential on-off model is approximated. 

The first approach that we shall consider uses a classical queueing 
theory model for the source. f Messages arriving at a terminal are 
viewed as customers arriving for service. An analogy is drawn between 
the length of the message (in bits or packets) and the service time of 



f The queueing theory model of a source in message switched networks is widely 
used, e.g., Itefs. 9 and 10. 
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a customer in queueing theory. Thus the exponential on-off source is 
approximated by messages with exponential length which arrive at 
a Poisson rate. 

The queueing theory model is well suited to sources that are not 
very active. (However the line into which these sources are multiplexed 
may be very active if many sources are connected to the ring.) When 
the source is more active the queueing model has some inherent in- 
accuracies. The queueing model implies that the distribution of time 
between messages is an exponentially distributed random variable 
and there is a nonzero probability of successive messages overlapping. 
In contrast for the exponential on-off model, the distribution between 
beginnings of successive messages is the convolution of two exponential 
distributions and the probability of message overlap is zero. 

The second approach to delay calculation uses a smoothed version 
of the traffic out of a source. Thus if the average number of bits per 
second out of a source is X bits/second, the calculation assumes that 
bits emanate from the source at a constant rate of X bits/second. This 
model of the source is meant to take up where the previous model 
leaves off, i.e., active sources. The effect of smoothing the bit flow will 
be less deleterious for more active sources. 

The analysis of message delay based on a modified M/G/l queue 1 " 
is based on work by B. Avi-Itzhak and P. Naor. 11 The line into which 
a message is multiplexed is viewed as a server that is subject to random 
breakdown. As the line is either idle or busy the server is operating 
or under repair. Four assumptions on probability distributions are 
necessary in order to carry out the analysis: messages arrive at a Poisson 
rate, the interval between message arrivals is independent of the 
message size, the duration of the line idle period is exponentially dis- 
tributed, and the lengths of line idle and busy periods are statistically 
independent. The distributions of the size of a message and of the 
duration of a line busy period are arbitrary. 

We take the arrival rate of messages from source i as X,- . The amount 
of time required to multiplex a message onto a completely free line 
at station i is denoted by the random variable S t . It is assumed that 
the duration of the idle period on the line is exponentially distributed 
with parameter A* [see equations (17) and (21)]. The random variable 
Li denotes the duration of the line busy period as seen by the source 
at station i. In Appendix B it is shown that the Laplace-Stieltjes 



t According to standard queueing theory notation an M/G/l queue is one where 
a single server accommodates customers that arrive at a Poisson rate with an arbi- 
trarily distributed service time per customer. 
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transform of T { , the delay suffered by a message at station i, is given 
by 

T ' {) ( ' 79,) X, -»- X..£ Wi (i»+ A? - A*£,,00) ' 

(24) 

where £ Li (v) and cC s ,(y) denote the L — S transforms of L, and £, 
respectively. By differentiating £ r ,(y) with respect to v and allowing 
v to approach zero, one obtains the following expression for the expected 
value of Ti, 

E[T { ] = E[S,} 0% 

, „,«., x,q + e*) a 
+ E[Si] 2[i - 7 ,^.(i + e,-)] 

+ £[L ' ] 2(1 + ei)[i -' 7,0,(1 + efl] ' (25) 

Now the mean number of packets per message is Q { . Since the line 
packet rate is C„ packets per second, the average time that it takes 
to multiplex a message onto a clear line is 

E[S,} = Qt/C, = J*- 

In the previous section we have made estimates of the line busy period 
which we have designated as 1/M* { . These quantities can be substituted 
into (25) yielding 

7,-e? , 7? 0-(i + e^'-'a + p at ) 



E[T { ] = 



where 



and 



M, M, 2[1 - 7,0,(1 + e*)] 



_1_ UW + g L< ) __ 

^ M* 2[i - 7.0,(1 + e*)] ' k ; 



Var (S.) 
P8i ~ E 2 (Sd 



Var (/,,■) 
PLi ~ E 2 (L.) 



The quantities /3. Sl . and f3 Li indicate the sensitivity of the calculation 
of delay to assumptions about probability distributions. In the next 
section sample calculations are presented in which we assume that 
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fj Li = p s . = l as would be the case if S { and L, were exponentially 
distributed. 

By taking a second derivative of equation (24) and letting v approach 
zero, the second moment of delay can be found. For brevity we have 
omitted this expression; however, calculations based on it will be 
presented in Section VII. 

£ Ti (v) can also be used to calculate the probability that a message 
has zero delay. Let us assume that lim u ^ m £ Li ( v ) — an d that S 
is exponentially distributed. Then 

Pr [zero message delay] = lim £ Ti (v) 



= (1 - Uf - ^d,) 



7M. 



(27) 



_7M. + A=L 

As mentioned previously, the second calculation of delay is predi- 
cated on the assumption that the data flow from the source is at a 
constant rate equal to the average rate from a source. Thus we assume 
that the source associated with station i generates data at a constant 
rate of r packets per second. Because the line is not continuously 
available to receive these packets buffering is required. It is the content 
of the buffer that is the key element in the calculation of delay. 

As stated earlier, the line traffic consists of alternate busy and idle 
periods. In the following it is assumed that the durations of the busy 
and idle periods are independent and exponentially distributed. The 
mean values of these quantities are known from the calculations in 
the previous section. Under these assumptions, constant input rate 
and exponentially distributed durations of line idle and busy periods, 
it can be shown 12 * that the probability density of the content of the 
buffer in packets of the ith station is given by 

/(£,) = KMB<) + (1 - KJai exp (-<*,£,), (28) 

where 5(-) is the Dirac delta function and K t and a, are related to 
the parameters of the system by 

K C„ - r ,-(l + e*) (99) 

Ki - (C„ - r,)(l + G*) ^ 



and 

Mf A* 



Ti C p — r, 



(29b) 



t In the reference, the solution was obtained for a hyper-exponential density. 
In the present work the density is exponential, which can be obtained from the 
hyper-exponen tial density. 
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Notice that K { is the probability that buffer i is empty. The average 
content of buffer /' in packets is 

E[Bi] = 1 ~ Ki ■ (30) 

The delay suffered by a packet is the amount of time that it spends 
in the buffer while waiting to be put on the line. From Little's Theorem 13 
the average delay of a packet is the average content of the buffer 
divided by the arrival rate: 

E[T t ] = EiBMu , i = 1, 2, ■ ■ • , N. (31) 

Substituting (29) and (30) into (31) yields 



^-Ut '■ 



Mf D - T.w.(l + 



(32) 



VI. RESULTS — AVERAGE DELAY 



This section is devoted to numerical examples of the foregoing 
results. We shall look at two contrasting modes of system operation, 
complete symmetry and complete asymmetry. In the symmetric case 
all stations transmit equally to all other stations and the destination 
matrix is given by equation (11). All stations in this case are precisely 
alike in their traffic characteristics. This symmetric traffic pattern 
may be encountered on a national ring which connects regional rings. 
The presumption here is that each of the regional rings to which it is 
connected receives and transmits the same volume of traffic. 

The asymmetric case models the situation where one station on the 
ring receives all of the output of the other stations. This singular 
station, in turn, distributes its traffic equally among the remaining 
stations. The asymmetric traffic pattern will be encountered in a ring 
which is composed of inquiry response users connected to a computer 
through a C-box. 

Let us begin with the symmetric ring. Calculations were made for 
two separate models of delay. Model 1 is based on the queueing model 
of delay that resulted in equation (26). The model of line busy period 
here is that which led to equations (16)-(19). The combination of these 
two models is best suited to the situation where there are many lightly 
loaded sources on the line which send short messages. Model 2 is a 
combination of the smoothed approximation to the source that led 
to equation (32) and the calculation of line busy period lengths that 
resulted in equation (23). 
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The estimates of average delay yielded by Models 1 and 2 are given 
in terms of average message lengths. Presumably, in a particular 
application, the duration of an average message is known and the 
delay in seconds due to line traffic can be evaluated. 

The results of the calculation for the symmetric ring are shown in 
Figs. 6-10. On Figs. 6 and 7 delay normalized to the average message 
duration is shown as a function of the number of stations for two 
different values of source utilization. The difference between these 
curves lies in the factor y which is equal to 1 on Fig. 6 and 1/30 on 
Fig. 7. The value of 1/30 for y corresponds roughly to sources with 
50 X 10 3 -bit-per-second active rate feeding into a 1.5 X lC-bit-per- 
second line. The results predict that, for y = 1/30, the ring can ac- 
commodate many stations with delay less than one average message 

length. 

Another view of system performance is shown on Figs. 8-10, where 
delay is shown as a function of source utilization for 10-, 50-, and 100- 
station rings when 7 = 1. We see from these curves that, for moderate 
line loading, delays are not large. For example, when the line utilization 
on a 50-station ring (see Fig. 9) is 0.5, the delay is less than two average 
message lengths. Also shown on Figs. 8-10 are the results of simulations 
which will be discussed presently. 
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Fig. 6— Normalized delay for u =0.1 and 0.01 (uniform inputs), y - 1. 
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Fig. 7 — Normalized delay for u = 0.1 and 0.01 (uniform inputs), y = 1/30. 
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Fig. 8 — Normalized delay for 10-station ring (uniform inputs). 
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Fig. 10 — Normalized delay for 100-station ring (uniform inputs). 
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A comparison of Figs. 10 and 1 1 shows the impact of different values 
of 7. On both figures delay is shown as a function of line utilization 
for 100-station rings. The difference between them is that y = 1 on 
Fig. 10 and 7 = 1/30 on Fig. 11. For U = 0.5 and 7 = 1 (Fig. 10), 
we find that the delay is two message units. In contrast, for 7 = 1/30 
and the same line loading (Fig. 11), the delay is near zero. When 7 = 
1/30, the line has very large capacity compared to the source data 
rates, even when the line is moderately loaded. 

It is convenient to use Model 2 to examine the asymmetric case 
of many users communicating with a single computer on the ring. 
Let u c denote the utilization of the computer and u, denote the utiliza- 
tion of the customer in each of the N stations connected to the computer. 
We assume that the computer distributes its traffic equally among all 
the other stations. Applying equations (7) and (9) the line utilization 
as seen by the ith station after the computer is 

Uf = y t u e S~£ + y,u s (i - 1), t - 1, 2, ... , N, (33) 

where y, and 7, are the ratios of source packet rate to line packet 
rate for the computer and the user, respectively [see equation (2)]. 
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Fig. 11 — Normalized delay for 100 and 500 stations (uniform inputs), 7 = 1/30. 
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Fig. 12 — Normalized delay at most critical station for N = 5 and 100 stations, 
7 = 1/30 (asymmetric case). 



In our computations we assume that the computar-to-user traffic is 
ten times the user-to-computer traffic (u c y c = 10 y,Nu„)S Equation 
(33) leads us to conclude that under this assumption the most critical 
station in the ring is the one right after the computer (i = 1) since 
the line traffic at this point is heaviest. 

Equation (33) may be used in conjunction with equation (32) to 
find the delay normalized to the line busy period at the most critical 
station. Because of the way the line busy period is calculated in Model 2, 
we may in turn normalize the line busy period to the activity time of 
the computer. The results are shown in Fig. 12 where delay, normalized 
to the computer busy period, is shown as a function of u, and u e and 
the number of stations as parameters. In Fig. 12 we have taken 7, = 
1/30 and 7r = 1. 



VII. SIMULATION RESULTS AND COMPARISON WITH THEORY 

An investigation of the system was carried out by means of simula- 
tion as well as analysis. Single rings comprising the one A-box and 
from 10 to 100 B-boxes were simulated. All of the simulation results 
were for symmetric rings in which each station has identical traffic 



t Studies of computer traffic support this assumption.' 
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characteristics and transmits equal portions of its traffic to the other 
stations. 

In the simulation, as in the analysis, our attention was focused on 
average delay. Nevertheless, as we shall see, the simulation yielded 
information on other characteristics of delay which can be compared 
to analytical results. 

The simulations attempted to mirror as closely as possible the actual 
operation of the ring. At each station a sequence of message and idle 
period lengths are chosen randomly from exponential distributions. 
As the messages arc generated they are assembled into an integral 
number of packets. Bit stuffing is used to round out the last packet 
in a message. The destination of each message is found by a random 
selection from N — 1 equally probable choices. 

After packets are generated they are given an initial time tag and 
sent to a buffer. The packets are stored in the buffer until they are 
multiplexed on the line. When a packet is multiplexed on the line, 
the time is noted and a difference in multiples of packets is taken with 
the initial time for each packet. It may happen, especially in a lightly 
loaded system, that a packet is multiplexed on the line immediately. 
In this case the time difference is zero. 

These time differences indicate the delay suffered by a packet due 
to line traffic. By noting the time difference for the last packet in a 
message, we have a measurement of message delay. Histograms of 
message delay were compiled and results drawn from these histograms 
will be presented in the sequel. 

A key step in the theoretical calculation was the estimation of line 
busy and idle period durations. Accordingly, in order to check the 
consistency of the models, measurements were made of line busy and 
idle period durations. These measurements will also be presented in 
the sequel. 

In order to keep the simulation effort within reasonable bounds, 
it was necessary to fix some of the parameters of the system. Thus for 
most of the data that follow the average message length was fixed 
at 100 bits. The messages were quantized into 125-bit packets (B„ = 
125). Header information was neglected (H = 0). Source utilization 
was varied by choosing appropriate idle durations. In this case 70 
percent of the messages are of one packet duration. As a check, selected 
simulations were run with different ratios of message length to packet 
size. 

As seen in the previous section, the theoretical models estimated 
average delays of one or two average message lengths for light to 
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moderate line loadings. For y = 1/30, estimates of average delay 
were small even for reasonably heavy line loadings. In Figs. 8, 9, and 
10 average delays yielded by simulation are shown in comparison 
with theoretical results. These averages were taken over all stations 
for each of the line loadings shown. In Fig. 8 we see that for the 10- 
station ring the theoretical estimates given by Model 2 are quite close 
to simulation results, differing by substantially less than a message 
duration. The estimates produced by Model 1 for moderate line loads 
are also fairly good, overestimating delay by about 0.3 message duration. 
Theoretical estimates compare well with simulation results for the 
50- and 100-station ring as well. In general both models overestimate 
the simulation delays somewhat. For line loadings below 0.5, Model l's 
estimates are within 0.5 message duration above simulation values, 
while Model 2's estimates are somewhat closer. For line loadings 
above 0.5, Model l's estimates are closer to simulation results. 

It is significant that delay estimates obtained for line loadings below 
the knee of the curve (U = 0.5 for the 10-station ring) are well within 
an average message duration. Since the system would probably be 
operated in this region, delay estimates for these line loadings are 
important. 

Each of the simulation points presented above is an average of 
between 7000 and 14,000 data points obtained through simulation. 
Along with average delay the standard deviation of delay was estimated. 
Estimates of the standard deviation of average delay were obtained 
by assuming that the standard deviation obtained from simulation 
was equal to the standard deviation of the underlying distribution 
of delay. The results showed that standard deviations of the averages 
are reasonably low. For example, for a 10-station ring at a line loading 
of 0.43, 10,000 data points were taken and the standard deviation 
was estimated to be 1.61 average message durations. The standard 
deviation of the mean is estimated to be 0.016 average message dura- 
tion which is less than 2 percent of the average delay. The situation 
is the same on the 100-station ring. At a line loading of 0.52 where 
7000 data points were taken, the standard deviation of the mean was 
estimated to be approximately 2 percent of the average delay. 

Average delay is clearly not the only characteristic of delay that is 
important in judging the performance of a system. Both the simulation 
and the analysis provided results on characteristics of delay beyond 
averages. As mentioned earlier, histograms of message delay were 
compiled by the simulation program. On Fig. 13 the cumulative prob- 
ability graph of message delay is shown. Also shown on Fig. 13 is the 
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Fig. 13 — Cumulative probability graph of queueing delay in a 100-station ring. 
Data compared with estimates from Model 2. 



probability distribution predicted by Model 2 from the same line 
loading and system configuration [see equation (28)]. The probability 
distribution for Model 1 was not readily obtainable. The cumulative 
probability graph shown here is typical of other system configurations 
and line loadings. It shows a nonzero component at zero delay with 
diminishing probabilities of larger packet delays. 

In order to summarize simulation results along these lines, three 
sets of numbers will be presented. They are probability of zero delay 
and a conditional mean and standard deviation of delay. The mean 
and standard deviation of message delay here are conditioned on the 
message having nonzero delay. 

In Table I are shown the probability of zero delay given by simula- 
tion and by Model 1 [equation (27)] and Model 2 [equation (29a)]. 
For the smallest utilization in the 10-station ring the error between 
simulation and theory is within 15 percent. For the other line loadings 
the estimate obtained from Model 2 is within 10 percent of simulation 
values while Model 1 considerably underestimates zero message delay 
probability. For the highest loading considered, U = 0.7 on the 100- 
station ring, neither model predicts the simulation value very well. 

Estimates of means and standard deviations of delay conditioned 
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Table I — Probability of Zero Message Delay 





17-Line 




Theory 








of Stations 


Utilization 


Simulation 


Model 1 


Model 2 


10 


0.10 


0.906 


0.833 


0.918 


10 


0.20 


0.80 


0.688 


0.833 


10 


0.43 


0.565 


0.418 


0.624 


10 


0.54 


0.464 


0.311 


0.514 


100 


0.235 


0.773 


0.622 


0.769 


100 


0.52 


0.454 


0.318 


0.4S5 


100 


0.54 


0.418 


0.264 


0.420 


100 


0.70 


0.252 


0.178 


0.304 



on nonzero delay for both simulation and theory are presented in 
Table II. Comparison of the theoretical estimates of the mean shows 
that for the lightest loading in both the 10- and the 100-station rings 
the analysis estimates the delay within 3 percent. For the heavier 
loadings, excluding the U = 0.7 loading on the 100-station ring, the 
error is about 20 percent of the simulated results for Model 1 and 
about 25 percent for Model 2. 

In summary, the data presented in Tables I and II indicate that a 
substantial portion of the messages (i.e., those that have no delay) 
are accounted for, and that the remaining messages have an average 
delay that is predictable, in the lighter loads, within 3 percent, For 
example, in the 10-station ring with a line loading of 0.1, 90.6 percent 



Table II — Mean and Standard Deviation for Delayed Messages 
(In Average Message Lengths) 





^-Line 


Simul 


ation 


Theory 






Standard 






A-Number 


Utiliza- 


Mean, 


Deviation 


Model 1 


Model 2t 


of Stations 


tion 


E[T | T > 0] 


a[T | T > 0] 


Mean Std Dev 


Mean 


10 


0.10 


1.18 


0.79 


1.19 1.22 


1.18 


10 


0.20 


1.34 


0.99 


1.44 1.46 


1.44 


10 


0.43 


2.06 


1.89 


2.33 2.48 


2.58 


10 


0.54 


2 . 66 


2.75 


3.12 3.34 


3.7S 


100 


0.235 


1.57 


1 .59 


1.61 1.61 


1.69 


100 


0.52 


2.95 


3.20 


3.13 3.16 


4.26 


100 


0.54 


3.29 


3 . 55 


3.78 3.81 


4.83 


100 


0.70 


6.48 


7.S9 


5.60 5.65 


10.72 



1 The density in this ease is exponential, and thus the mean and standard deviation 
are equal. 
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of the messages have no delay and the remaining 9.4 percent have an 
average delay of 1.18 average message durations. 

One observed difference between theory and simulation lies in the 
prediction of zero message delay by Model 1, which, for moderate 
line loadings, noticeably underestimates this quantity. There is evidence 
that the difficulty here may lie in the way that Model 1 treats quan- 
tization effects. For example, message lengths were taken to be ex- 
ponentially distributed when, in fact, they are geometrically distributed. 
As pointed out in the beginning of this section, most of the simulation 
results are carried out for an average message length of 100 bits and 
a packet length of 125 bits. Investigations of the effect of quantizing 
messages into packets are continuing. 

As mentioned earlier in this section, data were gathered on line 
busy and idle periods in order to check the consistency of our models. 
In Table III are shown the average durations of line busy periods for 
different loadings measured at one station in the ring. For comparison 
the average line busy periods predicted by Models 1 and 2 are also 
shown. Model 1 underestimates the length of the busy period by 
approximately 10 percent. The estimates given by Model 2 consistently 
overestimate the duration of the line busy period. For line loadings 
greater than 0.5, the estimate is poor. 

The average durations of line idle periods as found in simulation 
are shown in Table IV. Again Model 1 underestimates while Model 2 
overestimates. The error for Model 1 is higher than for the corresponding 
line busy periods. For line loadings greater than 0.5, Model 2 yields 
high estimates. 

If the data on the durations of line busy periods obtained from 



Table III — Line Busy Periods 1 " (Packets) 





[/-Line 


Simulation 


Theory 


No. of 




Standard 






Stations 


Utilization 


Mean 


Deviation 


Mean (Model 1) 


Mean (Model 2) 


10 


0.10 


1.58 


1.02 


1.53 


1.65 


10 


0.20 


1.81 


1.41 


1.70 


1.98 


10 


0.43 


2.66 


2.74 


2.24 


3.36 


10 


0.54 


3.20 


3.64 


2.80 


4.76 


100 


0.235 


1.94 


1.52 


1.78 


2.36 


100 


0.52 


2.96 


3.07 


2.80 


5.87 


100 


0.54 


3.32 


3.87 


3.24 


7.9 


100 


0.70 


5.08 


7.32 


4.49 


14.5 



+ The results shown here are averages for a single station. The line utilization is 
averaged over the entire ring. 
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Table IV — Line Idle Periods 1 " (Packets) 
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Utilization 


Mean 


Deviation 


Mean (Model 1 ) 


Mean (Model 2) 
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0.10 
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0.20 


7.22 


6. SO 


6. SO 


7.93 
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0.43 


3.40 


2.07 


2.S4 


4.48 


in 


0.54 


2.74 


2.24 


2.34 


3.95 


100 


0.235 


6.34 


5.71 


5.95 


7.70 


100 


0.52 


2.06 


2.75 


2.54 


5.85 


100 


0.54 


2.7S 


2.39 


2.3 


5.65 


100 


0.70 


2.05 


1.58 


1.95 


6.20 



t The results shown here are averages for a single station. The line utilization 
is averaged over the entire ring. 



simulation are used in the calculation of average delay the conclusions 
do not change substantially. Except for the U = 0.7 load point in 
the 100-station ring, the average delays predicted by Model 1 increase 
by less than 10 percent. Below line loadings of U = 0.5, the predictions 
of Model 2 decrease but are still fairly close to simulation points. 
Regions where the line loading is light to moderate are of greatest 
interest since it is most likely that systems would be operated in this 
region. For higher line loadings, delay in terms of average message 
lengths is large and small changes in loading lead to large changes 
in delay. 



VIII. CONCLUSION 

We conclude with a summary of our results. The analysis and the 
simulation of 10-, 50-, and 100-station rings show that for 7=1, i.e., 
source rate and line rate equal, message delay is less than 2 average 
message durations for line loading up to 0.5 of capacity. For line loadings 
greater than 0.5, delay increases substantially. For 7 = 1/30, i.e., 
a low-speed source feeding into a high-speed line, average delay in 
terms of average message durations is low for line loadings up to O.S. 
Analytical and simulation results indicate that there is a nonzero 
probability of zero message delay at all line loadings. The probability 
of zero message delay is greater than 0.5 for line loadings less than 0.5. 
Analysis and simulation indicate that the probability of a specific 
message delay decreases monotonically with the value of delay. The 
rate of decrease is similar to that of an exponential distribution (see 
Fig. 13). 
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For light to moderate line loads the theoretical models of average 
message delay show good agreement with the simulation results. For 
high line loads some of the approximations used in the theoretical 
models are weakened and agreement is not as good. However, even 
for heavy loading, there is meaningful agreement between analysis 
and simulation since both predict large delay and large shifts in delay 
with small shifts in line loading. 

Some of the discrepancy between analysis and simulation may be 
due to the fact that the analysis did not take into account quantizing 
effects. For example, line busy periods must be an integral number 
of packets in duration, whereas the analysis treated line busy periods 
as a continuous random variable. Investigation of these effects on 
message delay will be carried out in the future. 

The simulation program is capable of simulating systems where 
the traffic pattern is not symmetric and work will continue in this 
direction. 
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APPENDIX A 

Effect of Quantization 

In the text we have assumed that the duration of a source active 
period is exponentially distributed with parameter p ( (for source i). 
During the active period the source emits bits at a constant rate of 6, 
bits/second. The data is bundled into packets with B v bits/packet. 
The probability that there are j packets in a message is given by 

D t = / m, cxp (-M.0 dt, i = 1,2, ■•■ ,N 

= exp (-(; - 1)5 p m./&.-)[1 - exp (-B^/b,)]. 
Let 

A 4 cxp [-2WM; 
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then 

D, = A' _, (l - A). 
The mean number of packets per message is 

1/M< = £ jD, = z r-s— 7TT- 

frt 1 - exp (-BrUt/bi) 

APPENDIX B 

In an earlier paper Avi-Itzhak and Naor 11 found the average delay 
of a customer arriving at a server that is subject to random breakdown. 
Recently, using a similar line of reasoning, Avi-Itzhak 14 found the 
Laplace-Stieltjes transform of the density of this delay T\ to be: 



£ Ti ,(p) = (i - m - ye,) 

[Al|£ £ ..(tO - A? - v]£ Si (v + A? - AtJet.-fo)) 
\ t - v - \ t £ a ,(v + A* - A*i£ L ,(v)) 



(34) 



where £ Li (v) and «C Sl -( y ) & r e the Laplace-Stieltjes transforms of L< 
and «S, respectively. (Recall that in the main body of the text L, is 
defined as the duration of the line busy period at station i and S { is 
defined as the time required to multiplex a message on a free line.) 
In this derivation it is assumed that the line idle periods are expo- 
nentially distributed with parameter A*; . It is also assumed that 
messages arrive at a Poisson rate X, . 

In order for this result to be applicable to our problem, some modifica- 
tion of equation (34) is necessary. The expected value of T\ yielded 
by (34) when the line is entirely free for a long period of time is E[S t ]. 
But this is a delay that is due to multiplexing alone and is not a function 
of line traffic. Since we are interested in delay that is dependent on 
line congestion we remove this multiplexing delay. 

The delay described by equation (34) is the sum of two independent 
components, the waiting time and the residence time. The waiting 
time is the interval from when the message first arrives until it is first 
multiplexed on the line. The residence time is the interval in which 
the message is multiplexed on the line, including line busy periods 
during which the message is blocked. The L — S transform of the 
density of the residence time is given by 

£ Sl (v + At - A1£lM) 

= f° dt c~" t, [' dx (A '-ff - jt\{t - x)f 8i (x), (35) 

Jo *-0 ^O *>■ 
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where f ai (x) is the density of the message multiplexing time and ff t (x) 
is the fc-fold convolution of the line busy period. This expression is 
obtained by adding the message multiplexing times and all of the 
intervening line busy periods. Now if we simply add together only 
the line busy periods, removing the line multiplexing time, we have 

£ 5l .(A? - Al£ if (f)) 

= fW" t, f fc^7^«~" A, *ffI(0/*i(*)- ( 36 ) 

J i) k = u Jo '"'• 

Equation (34) becomes 

fi P ,(») - (1 - U* - yOt) 

\A*£ L ,(v) ~ A* - t'l&(Af - Al£,») (3?) 

X, - v - \<£ Si (v + Af - A ,-£ L ..(»)) 
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